Epidemic Diffusion Network of Spain: A Mobility Model to Characterize the Transmission Routes of Disease

Human mobility drives the geographical diffusion of infectious diseases at different scales, but few studies focus on mobility itself. Using publicly available data from Spain, we define a Mobility Matrix that captures constant flows between provinces by using a distance-like measure of effective distance to build a network model with the 52 provinces and 135 relevant edges. Madrid, Valladolid and Araba/Álaba are the most relevant nodes in terms of degree and strength. The shortest routes (most likely path between two points) between all provinces are calculated. A total of 7 mobility communities were found with a modularity of 63%, and a relationship was established with a cumulative incidence of COVID-19 in 14 days (CI14) during the study period. In conclusion, mobility patterns in Spain are governed by a small number of high-flow connections that remain constant in time and seem unaffected by seasonality or restrictions. Most of the travels happen within communities that do not completely represent political borders, and a wave-like spreading pattern with occasional long-distance jumps (small-world properties) can be identified. This information can be incorporated into preparedness and response plans targeting locations that are at risk of contagion preventively, underscoring the importance of coordination between administrations when addressing health emergencies.

Epidemic Diffusion Network of Spain: A Mobility Model to Characterize the Transmission Routes of Disease.

Human mobility drives the geographical diffusion of infectious diseases at different scales, but few studies focus on mobility itself. Using publicly available data from Spain, we define a Mobility Matrix that captures constant flows between provinces by using a distance-like measure of effective distance to build a network model with the 52 provinces and 135 relevant edges. Madrid, Valladolid and Araba/Álaba are the most relevant nodes in terms of degree and strength. The shortest routes (most likely path between two points) between all provinces are calculated. A total of 7 mobility communities were found with a modularity of 63%, and a relationship was established with a cumulative incidence of COVID-19 in 14 days (CI14) during the study period. In conclusion, mobility patterns in Spain are governed by a small number of high-flow connections that remain constant in time and seem unaffected by seasonality or restrictions. Most of the travels happen within communities that do not completely represent political borders, and a wave-like spreading pattern with occasional long-distance jumps (small-world properties) can be identified. This information can be incorporated into preparedness and response plans targeting locations that are at risk of contagion preventively, underscoring the importance of coordination between administrations when addressing health emergencies.

Identification of the first COVID-19 infections in the US using a retrospective analysis (REMEDID).

Accurate detection of early COVID-19 cases is crucial to reduce infections and deaths, however, it remains a challenge. Here, we used the results from a seroprevalence study in 50 US states to apply our Retrospective Methodology to Estimate Daily Infections from Deaths (REMEDID) with the aim of analyzing the initial spread of SARS-CoV-2 infections across the US. Our analysis revealed that the virus likely entered the country through California on December 28, 2019, which corresponds to 16 days prior to the officially recognized entry date established by the Centers of Disease Control and Prevention. Furthermore, the REMEDID algorithm provides evidence that SARS-CoV-2 entered, on average, a month earlier than previously reflected in official data for each US state. Collectively, our mathematical modeling provides more accurate estimates of the initial COVID-19 cases in the US, and has the ability to be extrapolated to other countries and used to retrospectively track the progress of the pandemic. The use of approaches such as REMEDID are highly recommended to better understand the early stages of an outbreak, which will enable health authorities to improve mitigation and preventive measures in the future.

Caveats on COVID-19 herd immunity threshold: the Spain case.

After a year of living with the COVID-19 pandemic and its associated consequences, hope looms on the horizon thanks to vaccines. The question is what percentage of the population needs to be immune to reach herd immunity, that is to avoid future outbreaks. The answer depends on the basic reproductive number, R0, a key epidemiological parameter measuring the transmission capacity of a disease. In addition to the virus itself, R0 also depends on the characteristics of the population and their environment. Additionally, the estimate of R0 depends on the methodology used, the accuracy of data and the generation time distribution. This study aims to reflect on the difficulties surrounding R0 estimation, and provides Spain with a threshold for herd immunity, for which we considered the different combinations of all the factors that affect the R0 of the Spanish population. Estimates of R0 range from 1.39 to 3.10 for the ancestral SARS-CoV-2 variant, with the largest differences produced by the method chosen to estimate R0. With these values, the herd immunity threshold (HIT) ranges from 28.1 to 67.7%, which would have made 70% a realistic upper bound for Spain. However, the imposition of the delta variant (B.1.617.2 lineage) in late summer 2021 may have expanded the range of R0 to 4.02-8.96 and pushed the upper bound of the HIT to 90%.

Assessing the effect of non-pharmaceutical interventions on COVID-19 transmission in Spain, 30 August 2020 to 31 January 2021.

BackgroundAfter a national lockdown during the first wave of the COVID-19 pandemic in Spain, regional governments implemented different non-pharmaceutical interventions (NPIs) during the second wave.AimTo analyse which implemented NPIs significantly impacted effective reproduction number (Rt) in seven Spanish provinces during 30 August 2020-31 January 2021.MethodsWe coded each NPI and levels of stringency with a 'severity index' (SI) and computed a global SI (mean of SIs per six included interventions). We performed a Bayesian change point analysis on the Rt curve of each province to identify possible associations with global SI variations. We fitted and compared several generalised additive models using multimodel inference, to quantify the statistical effect on Rt of the global SI (stringency) and the individual SIs (separate effect of NPIs).ResultsThe global SI had a significant lowering effect on the Rt (mean: 0.16 ± 0.05 units for full stringency). Mandatory closing times for non-essential businesses, limited gatherings, and restricted outdoors seating capacities (negative) as well as curfews (positive) were the only NPIs with a significant effect. Regional mobility restrictions and limited indoors seating capacity showed no effect. Our results were consistent with a 1- to 3-week-delayed Rt as a response variable.ConclusionWhile response measures implemented during the second COVID-19 wave contributed substantially to a decreased reproduction number, the effectiveness of measures varied considerably. Our findings should be considered for future interventions, as social and economic consequences could be minimised by considering only measures proven effective.

Perimeter Confinements of Basic Health Zones and COVID-19 Incidence in Madrid, Spain (preprint)

Background: A unique policy of perimeter closures of Basic Health Zones (small administrative health units) was implemented in the Autonomous Region of Madrid from September 21st 2020 to May 23rd 2021 to face the COVID-19 pandemic.Aim: To assess the impact of local perimeter confinements on the 14-days cumulative incidence of SARS-CoV-2 during the second wave of the pandemic in Madrid, Spain.Methods: We compare the errors in estimation of two families of mathematical models: ones that include the perimeter closures as explanatory covariables and ones that do not, in search of a significant improvement in estimation of one family over the other. We incorporate leave-one-out cross-validation and the choice of the best over 15 models at each step in our analysis for statistical signification.Results: The two families of models provided very similar estimations (correlation of the errors > 0.95 (±10-3 95% CI), difference in means of the errors < 1.2 (±0.7 95% CI) 14-days cumulative incidence), both for a 2 weeks and 3 weeks delay in observed cumulative incidence and also when restricting the analysis to only those Basic Health Zones that were subject to at least one closure during the time under study.Conclusion: Our analysis suggests that the perimeter closures by Basic Health Zone did not have a significant effect on the epidemic curve in Madrid, either 2 or 3 weeks after their activation.

Retrospective methodology to estimate daily infections from deaths (REMEDID) in COVID-19: the Spain case study.

The number of new daily infections is one of the main parameters to understand the dynamics of an epidemic. During the COVID-19 pandemic in 2020, however, such information has been underestimated. Here, we propose a retrospective methodology to estimate daily infections from daily deaths, because those are usually more accurately documented. Given the incubation period, the time from illness onset to death, and the case fatality ratio, the date of death can be estimated from the date of infection. We apply this idea conversely to estimate infections from deaths. This methodology is applied to Spain and its 19 administrative regions. Our results showed that probable daily infections during the first wave were between 35 and 42 times more than those officially documented on 14 March, when the national government decreed a national lockdown and 9 times more than those documented by the updated version of the official data. The national lockdown had a strong effect on the growth rate of virus transmission, which began to decrease immediately. Finally, the first inferred infection in Spain is about 43 days before the official data were available during the first wave. The current official data show delays of 15-30 days in the first infection relative to the inferred infections in 63% of the regions. In summary, we propose a methodology that allows reinterpretation of official daily infections, improving data accuracy in infection magnitude and dates because it assimilates valuable information from the National Seroprevalence Studies.